Renormalization group for phases with broken discrete symmetry near quantum critical points 
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We extend the Hertz-Millis theory of quantum phase transitions in itinerant electron systems to phases with 
broken discrete symmetry. Using a set of coupled flow equations derived within the functional renormaliza- 
tion group framework, we compute the second order phase transition line T c (6), with 6 a non-thermal control 
parameter, near a quantum critical point. We analyze the interplay and relative importance of quantum and 
classical fluctuations at different energy scales, and we compare the Ginzburg temperature To to the transition 
temperature T c , the latter being associated with a non-Gaussian fixed-point. 
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I. INTRODUCTION 

Quantum phase transitions in itinerant electron systems 
continue to ignite considerable interest. 1,2,3 These transitions 
are usually induced by varying a non-thermal control param- 
eter such as pressure or chemical potential at zero tempera- 
ture, so that no classical fluctuations occur and the transition 
is driven exclusively by quantum effects. In many physical 
situations, a line of finite temperature second order phase tran- 
sitions in the phase diagram terminates at a quantum critical 
point at T = 0. In such cases quantum fluctuations influ- 
ence the system also at finite temperatures altering measur- 
able quantities such as the shape of the phase boundary. Con- 
sequently, in a complete description of the system at finite T, 
quantum and thermal fluctuations have to be accounted for si- 
multaneously. 

The conventional renormalization group (RG) approach to 
quantum criticality in itinerant electron systems 4,5 relies on 
the assumption that it is sensible to integrate out fermionic 
degrees of freedom from the functional integral representa- 
tion of the partition function and then to expand the resulting 
effective action in powers of the order parameter alone. This 
approach has been questioned for magnetic phase transitions 
associated with spontaneous breaking of continuous spin ro- 
tation invariance, since integrating the fermions results in sin- 
gular interactions of the order parameter fieldi 2 ^ 3 - 

In this paper, we focus on quantum phase transitions to 
phases with broken discrete symmetry. We analyze quan- 
tum and classical fluctuations in the symmetry-broken phase 
with an Ising-like order parameter near a quantum critical 
point. Our calculations are based on a set of coupled flow 
equations obtained by approximating the exact flow equations 
of the one-particle irreducible version of the functional RG. 
Quantum and classical (thermal) fluctuations are treated on 
equal footing. The functional RG has been applied exten- 
sively to classical critical phenomena, 6 where it provides a 
unified description of 0(/V)-symmetric scalar models, includ- 
ing two-dimensional systems. The classical Ising universality 
class has been analyzed in Refs. 7,8,9. The functional RG 
is obviously also a suitable framework for treating quantum 
criticality^ In our approach, we can compute the RG flow 



in any region of the phase diagram, including the region gov- 
erned by non-Gaussian critical fluctuations. This allows com- 
parison between the true transition line and the Ginzburg line, 
which so far has been used as an estimate of the former. 3,5 
Specifically, we capture the strong-coupling behavior emer- 
gent in the vicinity of the transition line as well as the correct 
classical fixed-point for T c , including the anomalous dimen- 
sion of the order parameter field. 

Our results may be applied to commensurate charge- 
density-waves, Pomeranchuk transitions 11 ' 12,13 which sponta- 
neously break the discrete point group lattice symmetry, and 
to magnetic transitions with Ising symmetry, provided the 
transitions are continuous at low temperatures. 

The paper is structured as follows: in Section II we in- 
troduce Hertz's 4 effective action, adapted to the symmetry- 
broken phase, which serves as a starting point for the subse- 
quent analysis. In Section III we describe the functional RG 
method and its application in the present context, and sub- 
sequently derive the RG flow equations. In Section IV we 
present a solution of the equations in the case T — 0. Section 
V contains numerical results for the finite T phase diagram 
in the region with broken symmetry. Different cases are dis- 
cussed, distinguished by the dimensionality d and the dynam- 
ical exponent z. In particular, we compare the T c line with the 
Ginzburg line, thus providing an estimate of the critical region 
size. In Section VI we summarize and discuss the results. 



II. HERTZ ACTION 

The starting point of the standard RG approach to quantum 
critical phenomena in itinerant electron systems is the Hertz 
action. 4,5 It can be derived from a microscopic Hamiltonian by 
applying a Hubbard-Stratonovich transformation to the path- 
integral representation of the partition function and subse- 
quently integrating out the fermionic degrees of freedom. The 
resulting action is then expanded in powers of the order pa- 
rameter field, usually to quartic order. The validity of this 
expansion is dubious in several physically interesting cases, 
in particular for magnetic transitions with SU(2)-symmetry, 
since the integration over gapless fermionic modes can lead 
to singular effective interactions of the order parameter field, 
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which may invalidate the conventional power counting ! 2,3,14 
Such complications probably do not affect transitions in the 
charge channel and magnetic transitions with Ising symme- 
try. There are several indications that singularities cancel in 
that case, namely the cancellation of singularities in effective 
interactions upon symmetrization of fermion loops-l^ 6 - and 
the cancellation of non-analyticities in susceptibilitiesjiiiS We 
therefore rely on the usual expansion of the action to quartic 
order in the order parameter field and also on the conventional 
parametrization of momentum and frequency dependences, 
which leads to the Hertz action^ 5 - 



m 2 ^ J (2n) d 



+ U[<p] . (1) 



Here <p is the scalar order parameter field and <p p with p = 
(p, u>„) its momentum representation; a>„ = InriT with integer 
n denotes the (bosonic) Matsubara frequencies. Momentum 
and energy units are chosen such that the prefactors in front 
of jjpr and p 2 are equal to unity. The action is regularized 
in the ultraviolet by restricting momenta to |p| < Aq. The 
value of the dynamical exponent is restricted to z > 2. The 
cases z = 3, applicable to Pomeranchuk instabilities and Ising 
ferromagnetism, and z = 2, relevant for commensurate charge 
density waves and Ising antiferromagnetism, are of our main 
interest. Formally, the results of the paper may be applied to 
systems with arbitrary z > 2. A comprehensive discussion 
of the origin of the a>„ dependence of the action is given by 
Millis 5 - (see also Ref.[H). Eq. Q is valid in the limit ^ <K 1 
which is relevant here since the dominant fluctuations occur 
only in this regime. 

In the symmetric phase the potential U[<p] is minimal at 
= 0, and is usually parametrized by a positive quadratic 
and a positive quartic term. 4 ' 5 Since we approach the quantum 
critical point from the symmetry-broken region of the phase 
diagram, we assume a potential U[<f>] with a minimum at a 
non-zero order parameter <f>o: 
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where cf> and cf>' are functions of x and t with cf> = (/>(, + </>'. 
The parameter 5 - u <fr^/3 controls the distance from critical- 
ity. Approaching the phase boundary in the (S, T)-plane from 
the symmetry-broken phase gives rise to the three-point vertex 
V3 u 6, which generates an anomalous dimension of the order 
parameter field already at one-loop level. 

Although formally correct as a result of integrating out the 
fermions, the Hertz action (1) is not a good starting point for 
symmetry-broken phases with a fermionic gap, such as charge 
density wave phases or antiferromagnets. A fermionic gap 
leads to a suppression of the dynamical term (linear in fre- 
quency) in the action, since it suppresses low energy particle- 
hole excitations. If not treated by a suitable resummation in 
the beginning, this effect is hidden in high orders of pertur- 
bation theory. 20 We do not deal with this complication in the 



present paper. However, our results for the transition temper- 
ature should not be affected by a gap in the symmetry-broken 
phase, since it vanishes continuously at T c . 



III. METHOD 

To analyze the quantum field theory defined by S [</>] we 
compute the flow of the effective action F A [0] with approx- 
imate flow equations, which are derived from an exact func- 
tional RG flow equation 6, 21,22 ' 23,24 . The effective action T A [(f>] 
is the generating functional for one-particle irreducible vertex 
functions in presence of an infrared cutoff A. The latter is im- 
plemented by adding a regulator term of the form J ^<pR A (j) 
to the bare action. The effective action interpolates smoothly 
between the bare action S [</>] for large A and the full effective 
action T[(p] in the limit A — > (cutoff removed). Its flow is 
given by the exact functional equation 21 



d Ar 1 R A 



(3) 



where R A = d A R A , and F (2) [0] = 6 2 F A [(f>]/d(p 2 . In momen- 
tum representation (<p p ), the trace sums over momenta and fre- 
quencies: Tr = TYioj,, J For the regulator function /? A (p) 
we choose the optimized Litim cutoff 25 



tf A (p) = z(A 2 -p 2 )#(A 2 -p 2 ) , 



(4) 



where Z is a renormalization factor (see below). 

We approximate T A [(p] by a simplified ansatz with a local 
potential of the form Eq. (f2| characterized by cutoff dependent 
parameters u and <po (or, alternatively, 6), and a simple renor- 
malization of the momentum and frequency dependences of 
the quadratic part of the action via Z-factors, leading to an 
inverse propagator of the form 

G" 1 ^,^) = F (2) [0 = 0o] = Z w ^Lzp 2 +<5+/? A (p) , (5) 

IpK 2 

where the Z-factors depend on the scale A only. The regulator 
function R A (p) replaces Zp 2 with ZA 2 for |p| < A. 

Evaluating Eq. © for a momentum-independent field <p 
yields the flow of the effective potential U[<p] 



dfJJ[4>\ = - Tr — ;— : 

2 Z J^_ + Zp 2 + R A (p) + U"W 



(6) 



from which we derive the flows of the parameters <pu and u, 
following the procedure in Ref. 6. Viewing U as a function of 
p = j(f> 2 and using U'[po] = 0, we can write = j^U'lpo] = 
dfJJ'\po] + U"\po]d\po. Inserting dfJJ'\po\ as obtained by 
differentiating Eq. (|6]l with respect to p at p = po, and using 
U"[po] -\u, one obtains the flow equation for po 



OaPo = -Tr- — 

I \ 7 \u>„\ 



■Zp 2 +fi A (p) + | M p ] 



(7) 



The flow of u is obtained by differentiating Eq. (O twice with 
respect to p: 



3u 2 Tr ■ 



[Z w J#+Zp2 + flA(p)+§ M po] 



(8) 



Inserting the above flow equations into the A-derivative of 8 — 
tupo, we obtain the flow of 8. Below we will analyze the flow 
in terms of u and 8 instead of u and po. 

The Matsubara frequency sums in the flow equations © 
and (0 can be expressed in terms of polygamma functions 
*¥ n (z), defined recursively by M^+iCz) = T^,(z) for n 
0, 1,2,..., and T*o(z) = T'{z)IT{z), where T(z) is the gamma 
function. From the Weierstrass representation of the gamma 
function, r(z) _1 = ze yz + ^)e^ n , where y is the Euler 

constant, one can derive the relation 
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+ 2Wz+l) 



(9) 



Taking derivatives with respect to z yields expressions for 
sums involving higher negative powers of (\n\ + z) in terms 
of polygamma functions of higher order. 

For ^ A (p) we now insert the Litim function defined in 
Eq. ©. Following Ref. |g we neglect Z in ^ A (p), such that 
R A (p) = 2ZA©(A 2 - p 2 ) . The neglected terms would yield 
corrections to the flow equations that are linear in rj and there- 
fore irrelevant in the quantum part of the flow, where rj » 
(see below). In the classical part of the flow (for A — > 0) these 
corrections are small as compared to other included terms in- 
volving j] due to the presence of additional factors involving 
the interaction coupling and mass. Therefore they give rise to 
only minor corrections to the values of the fixed point coordi- 
nates and the values of the critical exponents. 

The (i-dimensional momentum integrals on the right hand 
side of the flow equations can be reduced to one-dimensional 
integrals, since the integrands depend only on the modulus of 
P 

Explicit dependences on A, Z-factors, and lengthy numeri- 
cal prefactors in the flow equations can be eliminated by using 
the following rescaled variables: 



p = 


Ipl/A, 


(10) 








T = 


ZA Z ' 


(11) 
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6 = 


ZA 2 " ' 


(12) 
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(13) 




FIG. 1 : Feynman diagrams representing the contributions to the flow 
equations dl4ll5U9t . 



The flow equations for 8 and S are then obtained as 

d8 



(flog A 



(1-2)8 



+ 4u 
+ \2u8 



d(l+8) 2 T 2 Jo 
d(l+dy T 3 Jo 



(14) 



du 



d\og A 



= (d - 4 + 2t]) u 



+ \2u 2 



1 1 



dd+Sf f 3 



-If 

f 3 Jo 



dp p d+3z - 7 V 2 (x) 



, (15) 



where x = 1 + T '(1 + 8)p z 2 , and 77 is the anomalous dimen- 
sion defined as 



1 = - 



dlogZ 
c/log A 



(16) 



To complete the system of flow equations one still needs to 
derive the evolution of Z and Z m , which parametrize the mo- 
mentum and frequency dependence of the propagator. Taking 
the second functional derivative of Eq. (|3), we obtain the flow 
equation for the propagator 

d K G- 1 (p) = 3u8 Tr [fl A (q) G 2 (q) G(q + p)\ , (17) 

where G(q) is given by Eq. (0. Here we skipped the contri- 
bution from the tadpole diagram, since it involves no depen- 
dence on momentum and frequency, and therefore does not 
contribute to the flow of the Z-factors. 

The momentum renormalization factor is given by 



Z= ^Ap[G- 1 (p,o*. = 0)] 



P=0 



(18) 



with A d = (2n)~ d S d .i, where S d -i = 2n d/2 /r(d/2) is the area 
of the (d- l)-dimensional unit sphere. 



where A p is the Laplace operator evaluated at constant cutoff 
function, that is, A p does not act on R A (p). Inserting Eqs. ( fTTl i 
and ( TT~8b into Eq. ( fl6b , and performing the frequency sum, one 



4 



obtains the flow equation for the anomalous dimension 



6 u5 r 
a (1 + 5) 5 L 



_8_i_££f- 5 f' 

' + 2^ 4 Jo 



dpp d+3z - n 



6p 4 (z - 2)(d + z- 4)f 2 <F 2 (;e) - 2/5 2+z (2(8 - d) 
(1 + 6 - p 2 ) + [(d - 14)(1 + 6) + 8/> 2 ]z + 3(1 + §)z 2 ) 
t^{x)-p 2z [2{p 2 -\) + ~6(z-2)+zf*¥ A {x) . (19) 



kept only the dominant term in the derivative expansion 6 and 
neglected the field dependence of the Z-factors. Furthermore, 
the simple parametrization of the effective potential Eq. (f2]l 
allowed us to substitute the partial differential equation © 
governing the flow of U[(f>] by the two ordinary differential 
equations ( fT4b and ( fT3T >. The latter approximation is equiv- 
alent to neglecting all higher order vertices generated during 
the flow. More sophisticated truncations have been applied in 
the context of the classical Ising universality class, where they 
lead to improved results for the critical exponents^ 



The contributions to the flow of 6, u and 77 are illustrated in 
terms of Feynman diagrams in Fig.Q] 

We now turn to the renormalization of Z„. We first consider 
the case z — 2, where Z„ can be related to the propagator by 



z " = ihr [ G ' (p = °' 2nT) ~ G ' (p = °' 0) J 



(20) 



Taking the logarithmic derivative with respect to A, inserting 
Eq. ( TP71 ). and performing the trace yields 



12 

t]o> = -r U $ T 
d 



T -{1+5) 



7X1 + c5) 3 



+ T-^ 2 1 + 



1 +6 



(21) 



From the numerical solution of the flow equation we observe 
that rju is small at all scales and vanishes for A — > 0. For 
example, in two dimensions 77^ varies between -0.033 and 
0.005 for u = 1 at T = e~ 4 and has practically no influence on 
the phase diagram. In three dimensions the values are at least 
one order of magnitude lower. For z > 2 one expects an even 
smaller r] w , since a larger z reduces the strength of fluctuations 
near the quantum critical point. Therefore we set 77^ = and 
Z w = 1 from now on. The scaling variable T then obeys the 
flow equation 



dT 



d\ogh 



in - z)t 



(22) 



In the flow equations (1 141 1 51 1 9b one identifies the classi- 
cal mean-field (involving only one power of 6 or 5), classical 
non-gaussian and quantum (involving f) terms. The quan- 
tum contributions vanish as the infrared cutoff tends to zero 
at constant non-zero temperature (T -1 <K 1). On the other 
hand the quantum terms dominate the high energy part of the 
flow, where T l » 1. Our framework allows for a continuous 
connection of these two regimes of the flow. The scale A c ; 
below which quantum fluctuations become irrelevant depends 
on temperature. It vanishes at T = 0. In the absence of a 
sizable anomalous dimension 77 in the quantum regime of the 
flow, one has 



A r l OC T 



Uz 



(23) 



as follows directly from the definition of f. It turns out that 77 
is indeed negligible down to the scale A c /, except in the case 
z = 2 in two dimensions (see Sec. 0. 

We emphasize that we have used a relatively simple 
parametrization of the effective action T A [(f)] . In particular, we 



IV. SOLUTION AT ZERO TEMPERATURE 

In this section we present a solution of the flow equations 
(1 1 41 1511 9b in the case of zero temperature, recovering earlier 
results for the critical exponents. We also provide an ana- 
lytic expression for the value of the control parameter 60 cor- 
responding to the quantum critical point. 

For this purpose we first take advantage of the asymptotic 
properties of the polygamma functions T'i(x) = x~ l + 0(x~ 2 ) 
for x » 1, T„ + i(x) = *¥' n (x), and evaluate the integrals in 
Eqs. (1 1 41 1511 9T > in the regime t <K 1. We introduce the vari- 
able 



- _ U_ _ A d 

V ~ T ~ 4^Z P Z W A 4 ~ d ~ z 



(24) 



and then pass to the limit T — > 0. Note that v does not depend 
on temperature. The resulting flow equations are 
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d\og A 

dv 
d\ogh 

'1 



(77-2)5 + 



I26v 



1 
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= (d + z - 4 + 2?7)v + 



(1 +6) 2 d + z-2 l+§d + z-2 
12v 2 1 



(1 + 6) 2 d + z - 2 ' 



Sv 



d(l+S) 4 



-a +5) 



(d-2)(z-2) 
3(d + z- 4) 



(1 +~5?~ 



d + z 



(25) 



The above expression for 77 is valid for d + z > 4. The flow 
equations for the case d — z — 2 are obtained by skipping 
the term involving (d - 2)(z - 2)/(d + z - 4) in Eq. ( 1251 
Non-gaussian contributions to the above equations originate 
from the quantum part of Eqs. (114115119b . The classical non- 
gaussian terms do not survive the limit T — > 0. 

Eqs. ( f25l ) have an infrared stable Gaussian fixed point in 
v = 0, 6 = 0, 77 = 0. Linearizing the flow equations around the 
fixed point, one obtains the solution 



6* + 



(d + z- 2) 2 
v(A) = (A/A ) 



((A/A ) d+Z - 2 - l) 



(A/A ) 



-2 



77(A) = 



(26) 



for d + z > 4. Here 6* and v* on the right hand sides denote 
the initial values of the parameters at A = Aq. In the marginal 
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case d + z — 4 one finds logarithmic convergence of v(A) and 
5(A) to zero. Expressing the order parameter 0o m terms of 
5* and v*, substituting the above solution, and taking the limit 
A — > yields 

0o <* V* 5 _ <*o . ( 27 ) 

where 

2A d A^' 2 ( ^ 
6 ° = — (d + z-2)i U (28) 

is the quantum critical point's coordinate. From Eq. (l27T i 
we read off the value of the exponent ft = I, consistent 
with mean-field theory. From Eq. (l26l > we can also straight- 
forwardly evaluate the correlation length using = 
limA^o S(A) = limA^o ZA 2 6, which yields 

€ = (S -SoT 1 ' 2 , (29) 

as expected within mean field theory. 

As anticipated, 4,5 the quantum phase transitions studied 
here behave similarly to classical transitions with effective 
dimensionality T> = d + z. For d > 2 and z > 2 one has 
T> > 4, leading to mean-field behavior governed by a Gaus- 
sian fixed point. 4 In the next section we turn to finite temper- 
atures, where non-Gaussian behavior occurs sufficiently close 
to the phase transition line. 

V. FINITE TEMPERATURES 

In this section we numerically solve the RG flow equations 
(1 141 1 51 1 91 > for finite temperatures. As already announced, the 
analysis is focused on the region of the phase diagram where 
symmetry-breaking occurs. In our notation this corresponds 
to sufficiently large values of the control parameter 6. First 
we treat the case z = 3 (in d — 2, 3) to which Eqs. (1 1 41 1511 9"T > 
are applied directly. The case z = 2, in which additional sim- 
plifications occur, is analyzed in Subsection B. In all numer- 
ical results we choose an initial cutoff Ao = 1, and an initial 
coupling constant u — 1 . 
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FIG. 2: (Color online) Anomalous dimension ?/(A) plotted as func- 
tion of s = — log(A/Ao) for different values of temperature in the 
case z = 3 and d = 2. The function rj(s) exhibits crossover from the 
mean-field value rj = 0, to the non-gaussian result j\ w 0.48. The 
crossover scale Ac is shifted towards smaller values of A (larger s) 
as temperature is reduced. 
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FIG. 3: (Color online) Quartic coupling f((A) plotted as function of 
s = - log(A/Ao) for different values of temperature in the case z = 3 
and d = 3. 



A. i = 3 

We first solve the coupled flow equations (1 141 1511 9b with 
the aim of determining the phase boundary T c (6), or, equiv- 
alently 6 C (T). To this end, for each given temperature we 
tune the initial value of 5 such that at the end of the flow (for 
A — > 0) one obtains the critical state with 6(A) — > 0. The 
tuned initial value is then identified as 6 C (T). Inverting the 
function 6 C (T) yields T c (5). In the variables 6, u this corre- 
sponds to seeking for such values of the initial 5, that both 
(5(A) and 5(A) reach a fixed point as the cutoff is removed. 

The flow of j] and u as a function of the logarithmic scale 
variable s — - log(A/Ao) is shown in two exemplary plots in 
Figs. [2] and [3] respectively. The flow is shown for various 
temperatures T, with 6 tuned to values very close to 6 C (T). 



The plateaus in Figs.[2]and[3]correspond to non-gaussian fixed 
point values of the flowing parameters. These values are not 
altered by the quantum terms in the flow equations. From 
Fig. |2j where d — 2, one reads off the value of the anoma- 
lous dimension rj * 0.48. The exact value from the Onsager 
solution to the Ising model is \. For the case d — 3 we find 
77 ~ 0.08 within our truncation, which is also about twice as 
large as accurate estimates of the exact value. 26 These results 
can also be compared to earlier classical 4 RG studies us- 
ing similar though more sophisticated truncations. In Ref.|7|a 
derivative expansion to order d 2 , retaining the general form of 
U[<p] and field dependence of wave function renormalization, 
was applied in dimensionality 1 < d < 4 yielding r\ m 0.08 in 
d — 3 and rj * 0.4 in d = 2. In Ref. 8 a derivative expansion 
to order <9 4 yielded rj as 0.033 in d — 3. To obtain scaling be- 
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havior in the range of a few orders of magnitude one needs to 
fine-tune the initial conditions with an accuracy of around 15 
digits. The breakdown of scaling behavior observed for very 
small A is due to insufficient accuracy of the initial value of 
6 and numerical errors. The plateaus are more extended as 
we go on fine-tuning the initial condition. Only exactly at the 
critical point true scale invariance manifested by plateaus of 
infinite size is expected. 

The figures also reveal the Ginzburg scale Ac at which non- 
gaussian fluctuations become dominant, such that the expo- 
nent i] attains a non-zero value. By fitting a power-law we 
observe Ac oc T c for d = 3 and Ac oc yfT~ c for d — 2 with non- 
universal proportionality factors. As expected, Ac vanishes at 
the quantum critical point, because at T — the effective di- 
mensionality D — d + z is above the upper critical dimension 
d c = 4. At finite temperatures, Ac is the scale at which u is 
promoted from initially small values (of order T) to values of 
order one. From the linearized flow equations one obtains 



Ac oc Tf 



(30) 



in agreement with the numerical results for d — 2 and d - 
3. Note that A G <s A r /, since A c / oc for z = 3, see 
Eq. d23l ). Hence anomalous scaling (finite rj) is indeed absent 
in the regime where quantum fluctuations contribute, and non- 
Gaussian fluctuations appear only in the classical regime. 

As already mentioned, the quantum contributions influence 
the flow only at relatively large A for T > 0. In particular, they 
do not alter any fixed-point values. However, at the beginning 
of the flow they dominate over the classical part and therefore 
are crucial for a correct computation of the initial value of 6 
leading to a scaling solution in the infrared limit. 

Results for the transition line T c (6) for d — 2 and d — 3 
are shown in Fig. [4] In both cases we recover the shape of the 
transition line as derived by Millis, 5 who used the Ginzburg 
temperature in the symmetric phase as an estimate for T c . 
Namely, we find 



for d - 2, 27 and 



(6 - 8 ) oc T c log T c 



(6 -6 ) 3 ' 4 



(31) 



(32) 



for d — 3. The exponent 3/4 in the three-dimensional 
case matches with the general formula for the shift exponent, 
ip = d+ l_ 2 for arbitrary z in dimensions d > 2. 5 Note that the 
phase boundary T c {5) approaches the quantum critical point 
with vanishing first derivative for d - 2 and with singular first 
derivative in the case d — 3. In Ref. 28 functional RG meth- 
ods in a different truncation were applied to calculate the shift 
exponent in the context of chiral symmetry breaking with the 
control parameter being the number of fermion flavours. 

An advantage of the present approach is that one can also 
follow the RG flow into the strong coupling regime, where 
non-gaussian critical behavior occurs. This in turn allows an 
estimate of the critical region's size as a function of tempera- 
ture or the control parameter 6. To evaluate the Ginzburg line 
in the symmetry-broken phase one solves the flow equations 
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FIG. 4: (Color online) The transition lines T C (S) obtained for z = 3, 
d = 2 (a) and z = 3, d = 3 (b). The phase with broken symmetry is 
located below the line. The phase boundary T c (6) obeys (6 — 5q) oc 
T c log T c for d = 2 and T c oc (S - So) 314 for d = 3 in agreement with 
the result by Millis. 5 The crosses indicate the Ginzburg temperature 
T G for various choices of 6. 



( 114I15I19D for fixed T and at different values of 5 > 6 C (T), 
observing the behavior of fixed-point values of the average 
order parameter <f>o (or, alternatively, the correlation length £) 
as 6 approaches 5 C . Typical results are plotted in Fig.|5]from 
which we read off the value of the exponent /3 describing the 
decay of the order parameter upon approaching the transition 
line 0q cc (S - 5cf. In the truly critical region (for 5- 5 C small 
enough) one obtains f3 ss 0.11 for d - 2 and ft ~ 0.30 for 
d = 3. These results come out close to the correct classical 
values 0.125 and 0.31, respectively. This is unlike the other 
critical exponents (77 and the correlation length exponent v) 
which within our truncation differ by factors close to 2 from 
their correct values. Indeed, as discussed in Refs. LZKa to ob- 
tain accurate values of the critical exponents in the Ising uni- 
versality class, and in particular in d = 2, one not only needs to 
consider the full partial differential equation governing the RG 
flow of the effective potential U[(f>], but also the field depen- 
dence of the wave function renormalization and higher orders 
in the derivative expansion of the effective action. 6 

From Fig. [5] we can extract the Ginzburg value 6c below 
which true critical behavior is found at the chosen tempera- 
ture T. Around #c, the exponent ft exhibits a crossover from 
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FIG. 5: (Color online) Order parameter cf> as a function of (S - 6 C ) at 
T = e~ 5 for z = 3, d = 2 (a), and z = 3, d = 3 (&). The exponent /3 
governing the decay of O upon approaching S c exhibits a crossover 
from a mean-field value /3 = 0.5 to a non-gaussian /3 « 0.11 for 
d = 2, and /3 ~ 0.30 for d = 3. The intersection of the straight lines 
determines Sa- 



lts mean-field value /3 = 0.5 to a non-gaussian value. In other 
words, 5c marks the boundary of the non-gaussian critical re- 
gion at a given T. At zero temperature, 5c coincides with the 
quantum phase transition point 6q, since there the fluctuations 
are effectively d+z > 4 dimensional, leading to mean-field be- 
havior. Several Ginzburg points in the 6—T plane are plotted 
as Tc{S) in Fig. |4] where they can be compared to the phase 
transition line. In three dimensions Tc and T c almost coin- 
cide, such that Tc provides an accurate estimate for T c . In 
two dimensions a sizable region between Tc and T c appears 
in the phase diagram. In that region, non-Gaussian classical 
fluctuations are present. 

We stress that accounting for the anomalous exponent 77 is 
necessary to describe the classical scaling regime, that is to 
obtain the plateaus in Figs. 2,3. Upon putting 77 = the scal- 
ing plateaus do not form. On the other hand, the shapes of the 
phase boundaries and the Ginzburg curves become very simi- 
lar at T — » in the present cases, where d + z > 4. Therefore 
our results for T C (S) are consistent with earlier studies, 5 where 
non-Gaussian critical fluctuations were not captured and the 
calculations were performed in the symmetric phase. 



In the case z — 2 the flow equations (1 141 1511 9b are signifi- 
cantly simplified, as all the integrals can be evaluated analyti- 
cally. One obtains 



dS 



dlogA 



= (77-2)5 + 2(5 



1 



+ MS 



d(\+ 6) 2 d 



d(l+S) 3 d 2W 



4 - , 



(33) 



da 



dlogA 



= (d - 4 + 277) u 



+ 3u 2 



4 1 4 - , 
d(l+6) 3 d 



(34) 



and for the anomalous dimension, 



u6 



d(l+6) : 



d 



d + 2 



¥400 



(35) 



where the argument of the polygamma functions is given by 
y = 1 + (1 + 6)/T . 

The procedure to evaluate the phase diagram and the 
Ginzburg line is the same as in the previously discussed case 
z = 3. In Fig.|6]we show results for the transition line T c {6) in 
two and three dimensions. We also show the Ginzburg tem- 
perature Tc for various choices of 5. In two dimensions, the 
transition line is consistent with the almost linear behavior 
(5 - So) T c log logTV/ log r c , derived previously for the 
Ginzburg temperature in the symmetric phase. 5 However, a 
sizable region with non-Gaussian fluctuations opens between 
T c and Tq. In three dimensions, T c {5) obeys the expected 5 
power law T c {5) oc(6— Sof with shift exponent if/ = 2/3, and 
the Ginzburg temperature is very close to T c for any 5. 



VI. SUMMARY 

We have analyzed classical (thermal) and quantum fluctu- 
ations in the symmetry-broken phase near a quantum phase 
transition in an itinerant electron system. The analysis is re- 
stricted to the case of discrete symmetry breaking, where no 
Goldstone modes appear. Following Hertz 4 and Millis^ we 
use an effective bosonic action for the order parameter fluctu- 
ations as a starting point. The renormalization of the Hertz ac- 
tion by fluctuations is obtained from a system of coupled flow 
equations, which are derived as an approximation to the exact 
flow equation for one-particle irreducible vertex functions in 
the functional RG framework. In addition to the renormaliza- 
tion of the effective mass and the four-point coupling, we also 
take the anomalous dimension 77 of the order parameter fields 
into account. In the symmetry-broken phase, contributions to 
j] appear already at one-loop level. Quantum and thermal fluc- 
tuations are captured on equal footing. The flow equations are 
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FIG. 6: (Color online) The transition lines T c (6) obtained for z = 2, 
d = 2 (a) and z = 2, d = 3 (b), The phase with broken symmetry 
is located below the line. For d = 2 the phase boundary r r (<5) is 
consistent with the relation (S - So) oc T c log log T c / log T c derived 
by Millis^ For d = 3 it follows the expected power law r c (<5) oc 
{S - So) 2,i . The crosses indicate the Ginzburg temperature Tq for 
various choices of 5. 



applicable also in the immediate vicinity of the transition line 
at finite temperature, where fluctuations deviate strongly from 
Gaussian behavior. 

We have computed the transition temperature T c as a func- 
tion of the control parameter 6 near the quantum critical point, 
approaching the transition line from the symmetry-broken 
phase. Explicit results were presented for dynamical expo- 
nents z — 2 and z — 3 in two and three dimensions. In all 
cases the functional form of T c {5) agrees with the behavior 
of the Ginzburg temperature above T c derived previously by 
Millis. 5 We have also computed the Ginzburg temperature Tq 
below T c , above which non-Gaussian fluctuations become im- 
portant. While Tc and T c almost coincide in three dimensions, 
a sizable region between Tq and T c opens in two dimensions. 

It will be interesting to extend the present approach to the 
case of continuous symmetry breaking. In that case Goldstone 
modes suppress the transition temperature more strongly, and 
more extended regions governed by non-Gaussian fluctuations 
in the phase diagram appear. 
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